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ABSTRACT 

We present the results of our multiwavelength study of two nearby Abell 
clusters that contain extended tailed radio sources. From our analysis of archival 
PSPC X-ray data, VLA radio observations, and new velocity measurements, 
we find strong evidence that these clusters are not relaxed. Significant X-ray 
substructure is revealed by examining deviations from circular models of the 
overall X-ray surface brightness maps. We report 130 new redshifts from the 
fields of A578 and A1569 and find that 41 galaxies belong to A578 and 54 belong 
to A1569. Moreover, we detect the presence of substructure in the velocity 
and spatial distribution of galaxies. We identify two interacting subclusters in 
A578 and two gravitationally unbound subclusters in A1569. Furthermore, we 
find that the radio sources in A578 are at rest with respect to their subclusters 
and that the Wide-Angled Tailed (WAT) source in A1569 has a small peculiar 
motion with respect to its subcluster (220 km/s). The standard model of ram 
pressure induced jet curvature requires relative velocities of the radio sources 
with respect to the ICM of ^ 1000 km/s. Therefore, the apparent peculiar 
velocities of the radio galaxies can not explain the bending of the jets/tails for 
the Narrow-Angled Tailed (NAT) radio source in A578 and for the WAT in 
A1569. We suggest that a cluster-sub cluster merger model provides a reasonable 
explanation for the X-ray and optical substructure detected in these systems. 
Numerical simulations of cluster-subcluster mergers predict a bulk flow of ICM 
gas that has sufficient dynamic pressure for bending and shaping the extended 
radio sources. Thus, a cluster-subcluster merger can provide the environment 
responsible for the observed cluster X-ray and optical morphologies and for the 
bending seen in the radio sources. 
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1. INTRODUCTION 

The study of substructure within galaxy clusters provides important constraints for 
cosmological models of the evolution of large-scale structures in the Universe. For instance, 
the amount of substructure detected in a cluster at a given epoch may be used to estimate 
the mean density of the Universe (e.g., |Richstone et al. 1992| ). Substructure is believed 



to be an indicator of an unrelaxed state in which the cluster is in the process of accreting 
mass via mergers or along large-scale filamentary structure. Observationally, substructure is 
detected in the spatial distribution of the cluster gas (i.e., in the X-ray surface brightness) 
or in the spatial and kinematical distribution of the galaxies. It has also been proposed that 
the presence of extended radio sources (i.e., Narrow-Angle tailed (NATs) and Wide- Angle 
tailed (WATs)) is another possible signature of a non- relaxed cluster ( Purns et al. 1994a . 



Burns et al. 19961) because close interaction between the radio emitting plasma and the 



turbulent ICM produced by a merger may be responsible for shaping the extended radio 
sources. In this paper, we test these hypotheses on two nearby Abell clusters that contain 
extended radio sources. 

Several methods have been developed for identifying and quantifying the amount of 
X-ray substructure present in galaxy clusters (i.e., Pavis fc Mushotzky 1993 . Mohr et al 



1993| , |White et al. 1994] ). In one of the pioneering studies of X-ray substructure, Jones 



& Forman (1984) found that up to 60% of galaxy clusters show evidence of substructure. 
Numerical simulations of cluster growth and evolution suggest that X-ray substructure 
originates from hot gas which is responding to a time-evolving gravitational potential 
( [Schindler fc Miiller 1993| , |Roettiger et al. 1996| , [Evrard 1990| , among others). The 



gravitational potential is continuously changing because dark matter and gas are being 
accreted through filaments and/or via mergers. Shocks, clumps, and elongations are 
typical observational characteristics of a cluster merger. Another effect of a merger is to 
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dramatically alter the gas temperature distribution as has been shown in recent ASCA 
and ROSAT observations (e.g., [Roettiger et al. 1995| , [Vlarkevitch et al. 1997|) . These 



results confirm the presence of shocks and temperature asymmetries as characteristics of 
non-relaxed clusters. 

Optically, a cluster that is still accreting matter and is not virialized may show 
strong signatures of substructure in the line-of-sight (LOS) velocity distribution and in the 
spatial (2D) distribution of its galaxies. Numerous statistical tests have been developed for 
detecting and quantifying optical substructure (see Pinkney et al. 1996 for an overview). 

Since the evolution of a cluster has significant effects on the intracluster medium 
(ICM) and on its galaxies, it is tempting to ask if these environmental changes affect 
the formation and development of radio galaxies. There is some evidence that low power 
extended radio sources (FR Is) are likely to be found in galaxy clusters that show evidence 
of substructure in the X-ray and in the optical ( Purns et al. 1994a| ). For instance, Pinkney 
et al. (1997) finds several examples of significant spatial-velocity correlations within WAT 
clusters associated with optically dominant galaxies at the optical/X-ray cluster centers. 
Gomez et al. (1997a) studied another sample of WAT clusters and found that 90% of their 
sample had significant X-ray substructure and X-ray extensions between the radio tails. 
They hypothesized that the X-ray substructure is evidence for recent mergers in these 
systems, which may explain the characteristic bending of the WAT radio jets. Finally, they 
proposed that ram pressure caused by the bulk motion flow of gas within a cluster-cluster 
merger is responsible for re-shaping the radio emitting plasma. 

Bliton et al. (1997) have recently found a possible correlation between the presence of 
NATs and X-ray substructure in a sample of optically similar Abell clusters with ROSAT 
PSPC data. They compared their results with a control sample of radio-quiet clusters and 
found that NATs are much more likely to be found in clusters with significant substructure. 
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They proposed a new model in which NATs are formed as a result of an increase in nearby 
ram pressure in cluster-cluster mergers. 

Finally, detailed studies of particular clusters show a very close relationship between 
the dynamical properties of galaxies, the local ICM, and the properties of extended radio 
sources. For instance, Rottgering et al. (1995) suggested that the radio sources in A2256 
(i.e., a radio halo and 5 NATs) were produced by the interaction of a disrupted ICM 
and the radio emitting plasma. They postulated that the ICM was stirred as a result of 
an on-going cluster-subcluster merger which was supported by strong optical and X-ray 
evidence. Moreover, Burns et al. (1995) explained the presence of an extended radio halo 
and NATs in A2255 as produced by a cluster merger. In their analysis, they compared the 
observed X-ray and kinematical properties of the cluster with hydro/N-body simulations of 
a cluster-subcluster merger. 

All of these results suggest that there is a correlation between the presence of extended 
radio sources and the dynamical state of galaxy clusters. We have begun to explore this 
intriguing possibility by performing a multiwavelength study of a sample of nearby Abell 
clusters. We are in the process of collecting redshifts and comparing our findings with X-ray 
images and numerical simulations of cluster evolution. In this paper, we report our results 
for two of the most interesting members of the sample: Abell 578 and Abell 1569. 

A578 has 2 prominent extended radio sources: a NAT and a double-lobed (Fat-Double) 
source (see Figure la). It is a relatively poor Abell cluster (richness class RC = 0) and has 
only 4 published velocities (Struble & Rood 1992 and a search through the NED archives). 
A1569 is also a RC=0 Abell cluster that contains a large WAT galaxy and an additional 
twin-jet source (see Figure lb). Pinkney (1995) has performed a kinematical study of 22 
galaxies with measured redshifts clustered around the WAT source in A1569 and found no 
significant evidence for non-normality. However, Burns et al. (1994a) and Gomez et al. 
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(1997a) found significant evidence for substructure from the X-ray surface brightness image. 
We present new optical and X-ray observations of these interesting clusters, and attempt to 
understand the relationship between the radio sources and the cluster environs. 

The paper is organized as follows. Section 2 describes the multiwavelength observations. 
In section 3, we discuss our optical and X-ray substructure analysis. In section 4, we 
analyze the dynamics of both clusters and compare our results with a cluster merger model. 
Finally, we summarize our conclusions in section 5. We use H G =75 km/s/Mpc and q o =0.5 
throughout. 

2. OBSERVATIONS 

2.1. Velocity Measurements 

Since we are interested in evaluating the dynamical state of these clusters, we need 
accurate positions and velocities for the cluster galaxies. Candidates galaxies were identified 
by eye from the POSS I plates and their coordinates measured from the CD-ROM Digital 
POSS I plates using the guide star plate solutions. The estimated relative accuracy of 
the digitized POSS I plates is ~ 0".5 . We found that the positions of the radio cores, as 
measured from the VLA radio maps, agree within 1-2 arcseconds. 

We identified a total of 178 candidate galaxies and 24 very bright stellar objects within 
~ one Abell radius (~ 2 Mpc or 20') in the direction of A578 and 70 candidate galaxies 
and 50 bright stellar objects within a similar region (~ 23'. 5) in the direction of A1569. 
We observed 123 objects in the direction of A578 and 66 candidate galaxies for A1569. In 
an effort to estimate the completeness and fairness of our sampling, we have compared 
the magnitude distribution of the observed objects with the overall distribution of all the 
candidate galaxies in the respective fields. We obtained rough magnitudes (with an error 
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~ 0.5) for these objects by using the CD-ROM Digital POSS I plates. We calibrated the 
candidate galaxies photometry by using the CCD Cousin R-band photometry of the radio 
galaxies in the clusters (Ledlow & Owen, 1995b). We find that the A1569 sampling is 90% 
complete down to m r =17.5. Our observations of A578 are 80% complete down to m r =16.5 
and 75% complete down to m r =17.0. 

All the spectra were obtained in 4 observing runs with the MX multifiber spectrometer 
mounted on the 2.3m Bok Telescope at the University of Arizona Steward Observatory. 
Table 1 shows a log of the observations. The MX multifiber spectrometer has 32, 2" fibers 
covering a field of view of 45'. There are 28 additional fibers used for obtaining sky spectra. 
The spectrometer has a 400 groove mm -1 grating which provides a dispersion of 2.8 A per 
pixel -1 and a spectral coverage from 3700 to 6950 A on an UV flooded Loral 1200x800 
thinned CCD. A more detailed description of the instrument can be found in Hill & Lesser 
(1988). 

We reduced the observations by following a very similar procedure to the one described 
by Hill & Oegerle (1993). We constructed a sky spectrum by combining ~ 28 sky 
spectra and then we subtracted this spectrum from each galaxy spectra. The scaling of 
the sky with respect to the galaxy spectra is based on the residual flux at the [OI] sky 
emission-line (A5577). Radial velocities were determined by using the IRAF task 'FXCOR' 
to cross-correlate the observed spectra with 20 templates. We used a similar criterion as 
the one described by Pinkney et al. (1994) to identify the best heliocentric velocity for the 
galaxies. 

Tables 2A and 2B contain the measured velocities for each cluster. The first column 
lists the galaxy ID number. Columns 2 and 3 list the epoch 2000 coordinates. Columns 4 
and 5 list the heliocentric velocities and errors. The errors were computed from a fit to the 
Tonry & Davis (1979) 'R' parameter (TDR) vs. velocity errors obtained from redundant 
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MX measurements of galaxy velocities ( |Pinkney et al. 1997 ). Column 6 shows the TDR 
values obtained from the cross-correlations. We have excluded spectra of very low S/N 
objects that had a low TDR number (< 3.5) or those that were clearly stars. We chose 
this TDR > 3.5 value cutoff because Pinkney et al. (1994) found that measurements 
with these TDR values would be generally good estimates of the redshift (~ 74% of 
correct cross-correlation peaks, Pinkney, private communication). Finally, the results of the 
cross-correlations were confirmed from visual inspection of the spectra. 



2.2. X-ray Observations 

The X-ray data for these clusters were obtained from the ROSAT PSPC pointed 
observations archive. A preliminary X-ray analysis of A1569 was reported by Gomez et al. 
1997a. We have performed a similar analysis for the A578 X-ray data. We used the Snowden 
et al. (1994) method because we are interested in the extended cluster X-ray emission. 
This method identifies four main contaminants within the X-ray data: Particle-induced 
Background (PB), After-pulse Background (AP), Scattered Solar X-ray Background (SB), 
and Long-term Enhancement Background (LTE). Next, it models each contaminant and 
corrects for them in each one of the Snowden ROSAT bands (R4, R5, R6, and R7) which 
roughly corresponds to the spectral range from 0.5 keV to 2.05 keV. Moreover, this method 
accurately creates exposure maps for every band. The final images consist of the summed 
maps over the range R4-R7. 



2.2.1. Temperature Fits 



The spectral analysis of the X-ray data was performed by using the XSPEO, 
'FTOOLS', and XSELECT' software in a similar fashion as in Gomez et al. (1997a). First, 
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we extracted two spectra from the data files with the same time filters as the ones used for 
the spatial analysis. One spectrum corresponds to the source region (with point sources 
removed) while the other spectrum corresponds to the background region. All the source 
regions had a radius of 500 kpc. Table 3 A shows the approximate center and radius of the 
source regions (500 kpc), and the inner radius and width of the background annuli. 

We rebinned the counts in the spectra such that every channel had at least 30 counts. 
Next, we fit the spectral profiles to Raymond-Smith ([Raymond fc Smith 1977| ) emission 
models (bremsstrahlung plus lines) in the energy range from 0.5 keV to ~ 2.05 keV. We 
assumed negligible internal absorption and a fixed cosmic abundance of 0.3 solar which 
is consistent with the values found by Edge & Stewart (1991a) for a sample of bright 
clusters observed by EXOSAT (their mean cosmic abundance is 0.32 with an rms of 0.18). 
We checked the quality of our 2 parameter fits (T and Njj) by comparing the resultant 
hydrogen column density with values obtained from a map of galactic HI (provided by F. J. 
Lockman, see |Stark et al. 1992]) . We found the fitted values to be in close agreement with 



those in the HI map. We then proceeded to fix N# (Table 3B, third column) to the Stark 
et al. (1992) values and fit only the temperatures with 'XSPEC 

In Table 3B, we report the results of the temperature fits. The fourth column lists the 
fitted temperature and the associated errors at the 68% confidence level for the \ 2 statistic. 
The fifth column lists the reduced x 2 from the fits and the number of degrees-of-freedom. 
The sixth column shows the temperatures from the C-statistic fits. These fits are based on 
the Maximum Likelihood method (Cash 1979); thus, they are considered more robust than 
the common \ 2 fits in cases in which the statistics are low for weak sources. 
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2.3. X-ray /Radio Overlays 

In order to compare the radio and X-ray morphologies, we first constructed final 
X-ray images in the range 0.5 keV to 2.05 keV. Then, we processed these images with an 
adaptive kernel smoothing algorithm ([Huang fc Sarazin 1997|) which readjusts the size of 
the smoothing length across the image to simulate a constant signal-to-noise ratio (~ 4) 
across the map. Next, we re-gridded and repositioned the X-ray images to match the VLA 
maps using the AIPS task HGEOM. These grids contained the same pixel size, field of 
view, and image center as the VLA maps. Grey-scale images of the radio emission were 
then overlaid onto the contours of the X-ray surface brightness. These X-ray/radio overlays 
are displayed in Figures la and lb. 



3. ANALYSES OF THE DATA 

3.1. Optical Analysis 

The analysis of the optical data consisted of the following steps. First, we determined 
the galaxy memberships in each cluster. Next, we performed ID, 2D, and 3D statistical 
tests in order to look for evidence of substructure. If we detected substructure within the 
cluster, we would attempt an objective partition of the galaxies and perform statistical tests 
on the resulting galaxy subgroups. 

In order to determine the cluster memberships, we first computed the biweight 
estimators for the scale (Sbi) and location (Cbi) ( Panese et al. 1980| , Peers et al. 1991]) 



of the velocity data. The main motivation for using these biweight estimators, which are 
analogous to the classical mean and dispersion, is that they are statistically robust. We 
computed the estimators errors by using the jacknife technique. Next, we applied the usual 
3<r-clipping criteria ( Yahil fc Vidal 1977Q for establishing the galaxy memberships in each 
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cluster. 

Once we were reasonably confident that we had eliminated most of the obvious 
interloper galaxies, we tested the data for substructure. For this task, we used Pinkney 
et al.'s (1996) BATCH program which performs a battery of ID, 2D, and 3D tests on the 
spatial and kinematical data. The ID section of the program includes the A, U, W, Bl, 
B2, B1-B2, I, DIP, KS, V, W 2 , U 2 , A 2 , TI, and Al-tests. Most of the 1-D tests compare the 
velocity distribution with a Gaussian distribution. Thus, the significance of the test reveals 
the probability that the two distributions are consistent; i.e. a low value for the probability 
means that the distribution is significantly non-Gaussian. Other tests analyze the shape of 
the distribution (i.e., TI test). The software also provides an average of the probabilities for 
non-normality, skewness, and kurtosis. 

The 2D tests are the Angular Separation Test (AST), the Fourier Elongation Test 
(FEL), the j3 test, and the 2D Lee statistic. The 2D Lee statistic uses reshuffling of 
the galaxy positions to estimate the significance of any substructure. Therefore, a low 
probability indicates significant substructure. 

Finally, the program uses the following 3D tests: Lee 3D, A (Dressler-Shectman or DS), 
e, the a tests (West-Bothum or WB), and the a variation test (WBv). The significance of 
these tests is determined by using a Monte-Carlo method. Thus, they test how many times 
the observed measurement is greater than the simulations based on the reshuffling of the 
observed data. This implies that low probability numbers identify significant substructure. 



Whenever we detected statistically significant substructure, we used the KMM ( [Bird 



et al. 1992 , pird 1993| , Pavis et al. 1995p algorithm to objectively partition the data 



in sub-groups. Briefly, this algorithm models the velocity distribution and the projected 
spatial distribution of the galaxies as composed by a user-defined number of sub-groups. 
This is the so called mixture modelling technique. Each one of the sub-groups is in turn 
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modelled as a Gaussian. The program finds the best fit by assigning different galaxies to 
different groups until the likelihood statistic is maximized. This program needs a first 
guess for the centroid (3D) position of the subgroups, their relative sizes (3D), and their 
estimated errors. As an additional check, the method also provides an estimate of the 
formal probability of rejecting the hypothesis that one Gaussian would fit the data. 



3.1.1. A578 Galaxy Velocities 

In Figure 2a we show a histogram of the velocity distribution for the galaxies measured 
in this field. The main cluster of 41 observed galaxies (which includes the two radio galaxies 
of interest) is located at ~ 27,000 km s" 1 while a background group of 17 galaxies is located 
at ~ 40,000 km s -1 . We have identified the NAT radio galaxy 0720+670 with our galaxy 
#5 and 0719+670 with galaxy #1 in Table 2A QOwen fc Ledlow 1997Q . The values for the 



biweight estimators and other data for A578 are listed in Table 4. 

Figure 2b shows a velocity-coded 2D plot of the spatial distribution of the galaxies in 
this field. The background galaxies appear to be uniformly distributed spatially whereas 
the galaxies clustered at 27,000 km s _1 are more centrally concentrated. 

Figure 2c shows a histogram of the velocity distribution of the accepted galaxies for 
A578 (these are all the observed galaxies clustered at 27,000 km s -1 that survived the 
3a clipping). This distribution does not appear to be Gaussian. In order to quantify the 
non-normal character of this distribution and to detect other signs of substructure, we ran 
the statistical tests described by Pinkney et al. (1996). Table 5 shows the results of the 
substructure analysis. The ID tests indicate possibly significant kurtosis (average = 9%) 
and a significant asymmetry (AI=5%). Although there is no evidence for substructure in 
the 2-D spatial distribution of the galaxies, we find evidence for significant 3D substructure 
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(2%) from the a test. Figure 2d shows the existence of an apparent spatial segregation 
between the low velocity galaxies and the high velocity galaxies. One group of galaxies 
is concentrated in the South with velocities < 27,000 km/s while the other galaxies with 
velocities > 27,000 km/s are predominantly located in the North. It is surprising, however, 
that the other 3D tests do not find any significant evidence of substructure. However, 
Pinkney et al. (1996) pointed out that these 3-D tests have widely varying sensitivity to 
different substructure dependent on the viewing angle relative to the merger axis and the 
epoch of the merger. 

Because we found significant ID substructure and possible 3D substructure, we used 
the mixture modelling technique in order to identify the subgroups present in the data. 
Table 6A shows the results of applying the KMM algorithm to A578. The method splits 
the galaxy population into two groups that we labeled North (18 members) and South (23 
members) due to their centroid positions. Most of the galaxy assignments which are listed 
in Table 6B had a formal probability > 95% to belong to their listed subgroups. Only 
galaxies #15, #20, and #27 had marginal probabilities > 80%. The observed distribution 
is consistent with a single Gaussian distribution at the very small 0.008% level. This 
significance might not mean that each galaxy assignment is correct per se, but it justifies 
our effort to characterize the components responsible for the substructure. To verify that 
the results were reliable, we performed the same analysis by changing the initial guesses 
for the sub-groups centroids and errors as suggested by Davis et al. (1995). The algorithm 
converged to the same results regardless of the guess for the initial parameters. Figure 2e 
shows a sub-group coded 2D spatial distribution of the galaxies from the KMM partitioning. 
It is reassuring to note that the new assignments show a distribution that is very similar to 
the naive velocity-based splitting depicted in Figure 2d. 

Figure 2f shows two contour plots overlaid onto a grey-scale image of the cluster X-ray 
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emission. The contour plots represent the galaxy surface density distribution of the two 
subgroups identified by the KMM algorithm. Both density distributions were smoothed 
with an adaptive kernel algorithm described by Beers et al. (1991). We see several 
interesting correlations between the X-ray and optical data. First, both galaxy subgroups 
appear to be elongated. The southern subgroup shows a E-W elongation while the northern 
subgroup exhibits a N-S elongation. Second, these elongations in the galaxy distribution 
appear to trace the cluster X-ray emitting gas. Finally, the centroid of the southern group 
is centered on a X-ray peak while one of the secondary peaks of the northern group is 
located close (1'.2 or 110 kpc) to a local peak of X-ray emission. 

Where are the radio galaxies located? Galaxy #1 is located at the center of the 
southern group and the NAT is located close to one of the subpeaks of the southern 
group (compare Figures la and 2f). Kinematically, they appear to be at rest within their 
subgroups. We estimated the significance of the radio galaxy velocities using the analysis of 
Teague et al (1990). They introduced the Sy statistic defined as follows: 



Sy = lCBI ~ V [ (1) 

where v is the galaxy velocity and <7i and <r 2 correspond to the errors in the cluster location 
and galaxy velocity. A value of Sy < 1 would mean that the galaxy is at rest and Sy > 2.0 
corresponds to a marginally significant peculiar motion. When we analyze the kinematics 
of the radio sources with respect to the overall cluster, we find that galaxy #1 (v pec =622 
km/s and S v — 4.0) and galaxy #5 (v pec =935 km/s and S v = 6.0) have significant peculiar 
velocities. However, after we correct for the substructure present in the systems, we find 
that these radio galaxies are in fact at rest with respect to their subgroups (i.e., Sy < 
1 when we only consider the velocities within each subgroup). This poses an interesting 
question: how is the southern tailed source bent if the radio galaxy is at rest? 
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3.1.2. A 1569 Galaxy Velocities 

Figure 3a shows a histogram of the velocity distribution for the galaxies measured in 
A1569. Almost all of the observed galaxies are clustered at ~ 22,000 km s _1 . There are 
only 2 galaxies located outside this velocity range (~ 31,500 km s -1 ). We have identified 
the two radio sources associated with A1569 with the galaxies clustered at ~ 22,000 km 
s _1 . Galaxy $4 corresponds to 1233+169 and galaxy #2 corresponds to the WAT radio 
source 1233+168 ( pwen fc Ledlow 1997] ). The values for the biweight estimators are listed 



in Table 4. Note the unusually high scale value (~ 1618 km/s) for a richness class cluster. 

Figure 3b shows a histogram of the velocity distribution of galaxies centered at 22,000 
km s _1 . This distribution is clearly non-Gaussian and shows the presence of at least two 
groups of galaxies. One relatively compact group appears to be clustered at ~ 20,700 km 
s _1 and the other more dispersed group is concentrated at ~ 24,000 km s -1 . This explains 
the very large value of the biweight scale for this cluster. However, are these two groups 
also separated spatially? In Figure 3c we show a velocity-coded 2D spatial distribution 
of the galaxies which indicates that the two groups are indeed spatially segregated. The 
low velocity galaxies are concentrated in the S-W area near 1233+168, whereas the high 
velocity galaxies are concentrated towards the N-E area near 1233+169. 

As expected, the statistical analysis finds strong evidence for substructure (Table 7). 
There is significant kurtosis (< 0.1%) and non-normality (3%). Moreover, the TI and AI 
reflect the velocity contamination present due to the N-E group. There is also evidence 
for 2D substructure. The /3 test reveals the asymmetry in the distribution of galaxies and 
the Lee 2D test points towards the apparent bimodality of the spatial distribution of the 
galaxies. Furthermore, the Lee 3D and e tests also show evidence for substructure. 

The KMM algorithm successfully partitions the galaxies into two groups. The basic 
data for the subgroups and the galaxy assignments are listed in Table 8A and Table 8B. The 
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mixture modelling technique assigns 25 objects to the N-E group and 29 objects to the S-W 
group (Figure 3d). Galaxy #51 had the lowest probability of group membership (99%). 
Moreover, the method finds that the distribution is consistent with a single Gaussian at 
< 0.1% level. Kinematically, the KMM method also separates the galaxies based on their 
velocities (Figure 3e). The galaxies located at the N-E have the highest velocities whereas 
the galaxies at the S-W have the lowest velocities. 

Figure 3f shows an overlay of the X-ray emission from the cluster (gray scale) onto a 
contour plot of adaptively smoothed galaxy surface density. As for the case of A578, we 
find that the galaxy density peaks are coincident with X-ray clumps and that both maps 
show similar elongations. Interestingly, the two subgroups appear to be connected by gas 
and by a tenuous trail of galaxies. 

We find that the WAT galaxy has a marginally significant (<Sy=2.2) peculiar motion 
(215 ± 105 km/s) with respect to the dynamical center of its subgroup. Note that this 
velocity, however, is insufficient to account for the bending of the jets. 

3.2. X-ray Analysis 

Another probe of the dynamical state of the cluster is the X-ray emission from its 
intracluster gas. Any substructure in the gas is probably caused by an evolving gravitational 
potential or deviations from hydrostatic equilibrium. Moreover, a study of the physical 
properties of the gas will allow us to determine if the cluster gas and the radio emitting 
plasma are in pressure equilibrium. 
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3.2.1. X-ray Substructure 

Figures la and lb reveal a very interesting spatial correlation between the position 
of the extended radio sources and the presence of X-ray clumps for our two clusters. In 
these Figures we note that all the radio sources appear to coincide spatially with clumps of 
X-ray emission. This effect was first noted by Burns et al. (1994a) from lower resolution 
Einstein images. In order to quantify the significance of these X-ray clumps, we performed 
a substructure analysis of the X-ray emission from these clusters. 

Our analysis for A1569 has been previously presented in Gomez et al. (1997a) and is 
similar to the method used by Davis & Mushotzky (1993). In short, our method identifies 
substructure by comparing the actual X-ray emission from the cluster with a circular model. 
Next, we estimate the significance of the substructure by using Poisson statistics. We have 
removed any obvious point sources during our analysis. Therefore, all the structure that we 
list in Table 9 as significant is not produced by point sources. 



3.2.2. Intracluster Gas Properties 

In order to characterize the general morphology of the clusters, we fitted f3 models 
(e.g., [Sarazin 1986| ) to the background subtracted profiles. The profiles were generated by 



azimuthally averaging concentric annuli with a bin size = 30" using Poisson statistics for 
the errors (background + sources counts added in quadrature). We used circular annuli in 
order to create the profiles and masked all the point sources in the region of interest. The 
best fit values for A578 are (3= 0.6 ± 0.1, a core radius of 120 ± 8 kpc, with a reduced \ 2= 
0.6. For A1569, Gomez et al (1997a) report values of (3= 0.6, a core radius of 182 ±110 
kpc, and a reduced x 2= 0.5. Although the quality of the fits are marginal, these fits are a 
convenient parameterization of the cluster properties. 
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We used the deprojection technique of Arnaud (1988) to estimate the density, 
temperature, and pressure profiles for the cluster gas. This method assumes a Raymond- 
Smith (1977) thermal spectrum, that the gas is in quasi-hydrostatic equilibrium, and 
requires an initial estimate for the temperature in the outermost boundary of the profile 
and the scale of the X-ray emission from the f3 model fit (/3 and r c ). 

The deprojection technique also needs an estimate for the velocity dispersion. In our 
case, we replaced the velocity dispersion by the biweight estimator of the scale. We used 
380 km/s for A578 (southern clump of galaxies) and 433 km/s for A1569 corresponding to 
the southern X-ray clump. The densities determined from the deprojections, along with the 
measured temperatures, were used to estimate the ICM pressures listed in Table 10. 

Finally, we compared the thermal pressures from the deprojection to the minimum 
pressure in the radio tails (e. g., P'Dea fc (J wen 1985| ). We used the radio images from 



the VLA 20-cm radio survey of Abell clusters ( Pwen fc Ledlow 1997| ) for this analysis. 
The minimum internal pressures were measured at the middle of the diffuse portions of the 
radio tails. For our calculations, we assumed a filling factor of 1, a power-law spectrum 
with a = 1 (F u oc u~ a ) with lower and upper cutoffs of 10 MHz and 100 GHz, a transverse 
magnetic field {sinO = 1), cylindrical symmetry, and an approximate equipartition of energy 
between the relativistic electrons and protons. 

Table 10 lists the results for each cluster. The tails of the two radio sources in A578 
appear to be in rough pressure equilibrium (within a factor of 2). However, the thermal 
and radio plasma pressures for A1569 are quite different. The tails of the southern WAT 
appear to be underpressured by more than a factor of 10 with respect to the surrounding 
ICM. A similar trend has been found in other cluster radio sources such as A2255 ( [Burns et| 



al. 1995| ) , Hydra A Qlaylor et al 1990]) , and a sample of cluster sources ( |Feretti et al. 1990| ) . 



The presence of thermal gas accreted from the ICM in the tails (i.e., entrainment) has 
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been suggested as an explanation for this effect because the entrained intracluster gas will 
contribute its own thermal pressure inside the radio tails which is not taken into account 
by the minimum pressure computation described above. 



3.2.3. Jet Bending 

In this section, we discuss a feasible mechanism for the bending of the jets based on the 
O'Dea (1985) and Doe et al. (1995) analysis of cluster NAT radio sources. In their studies, 
they modeled the NATs as 'cold' jets (Mach numbers 3> 1) that are bent due to ram 
pressure from interaction with the ICM. Thus, they assume that the jet can be modeled as 
a laminar relativistic fluid that obeys the relativistic and time independent Euler's equation 



Pi V h 2 „ PlCMVrel 2 ^ 



R h 



i 



where Vj is the jet velocity, pj is the jet density, f3j = Vj/c, 7 = (1 — /3 2 ) -1 / 2 , R is the radius 
of curvature of the jet, v re i is the relative velocity of the radio source with respect to the 
ICM,/i is the scale length over which the ram pressure is effective, and picu is the gas 
density. 

Next, they assumed that the initial jet kinetic energy is converted to the observed 
radio luminosity with an efficiency of e. This yields 



Vrel 



L rad f3j-fh ( r 



^)V3, (3) 



^ (7 - l)picM^rj 2 cR r if ' 

where r,-, and r^j are the average, initial, and final jet radius and L rad is the radio 
luminosity. As mentioned by Doe et al. (1995), this expression is very useful because it 
involves quantities that are easily determined from the multiwavelength observations of the 
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clusters (rj, Picmi L ra d). However, we still have to make some assumptions for e and (5j. 

As a first step in our analysis, we calculated the radio luminosities for the jets in A578's 
NAT and the Al569's WAT. These luminosities were computed from the flux densities by 
assuming a spectral index of 0.5 for the jets and a bandwidth from 10 MHz to 15 GHz as in 
O'Dea and Owen (1987). We obtained L A578 =0.83 x 10 41 erg/s and L A1569 =4.8 x 10 41 erg/s. 

Figure 4 shows plots of v re i and jet number density as a function of (3j for the southern 
A578 NAT radio source and for the southern A1569 WAT radio source. This Figure 
illustrates that the jet is more difficult to bend at high densities and/or if the jet is slow. 
On the other hand, a fast and/or tenuous jet is easier to bend. 

In order to further constrain the different parameters governing the jet bending, we 



have examined numerical simulations of the disruption of jets ( [Zhao et al. 1992]) ; these 
simulations show that a jet would be disrupted at a characteristic jet disruption length 
which is a function of the jet Mach number. By examining the radio data from these 
sources, we have determined that both jets are disrupted very close to their host galaxy. 
Thus, it is safe to assume that the jets propagate without disruption from some 20-30rj 
(these are their jet disruption lengths). This will imply Mach numbers between 4 and 5 (see 
Figure 7 in Zhao et al. 1992). Moreover, if we assume that the typical sound speed of the 
ICM surrounding the jets is ~ 500 km/s, we find j3j values in the range ~ 0.007 — 0.008. 
These (3j values constrain the possible relative velocities (Figure 4) between 400 - 2500 
km/s. This is the range of relative velocities of radio sources with respect to the ICM that is 
responsible for the observed radio morphology. But, how could these high relative velocities 
arise if the radio galaxies are at rest with respect to the other galaxies in their subgroups 
(see section 3)?. We will present a possible solution to this puzzling problem in the next 
section. 
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4. DISCUSSION 

Our analysis of the X-ray and optical data confirms our suspicions that these clusters 
are not relaxed systems. It is important to note that we initially chose these clusters only 
because they contain extended radio sources. In what follows, we describe our general 
impressions for each cluster and discuss models for the jet bending. 

4.1. A578 

From our multiwavelength observations of the A578, we find significant substructure in 
the gas density, as revealed by the X-ray surface brightness maps, and in the distribution 
of the cluster galaxies. When we compare the results of our X-ray and optical analysis, we 
discovered very interesting properties. First, the velocity dispersion of the whole system 
(793 km/s, Table 4) does not agree with the overall X-ray temperature of the gas. We 
estimate that a velocity dispersion of 500 km/s would be in better agreement with the X-ray 
temperature of this cluster according to the T-cr relationship for clusters of galaxies ([Edge 
fc Stewart 1991 b|) . Interestingly, the optical subgroups identified by the KMM algorithm 
have velocity dispersions (~ 400-500 km/s) that are in much better agreement with the 
cluster X-ray temperature. Second, these galaxy subgroups appear to coincide spatially 
with X-ray clumps (Figure 2f). 

We believe that these correlations and the other cluster properties that we have 
described in previous sections indicate that this cluster has undergone a recent cluster- 
subcluster merger. This merger will not only explain several of the X-ray and optical 
properties in this cluster, but it will also offer a solution to the puzzling problem of the 
NAT bending. 

In order to explore this merger model further, we have performed numerical simulations 
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of cluster-subcluster mergers and compared them with our observations. The simulations 
were generated with a code that is a hybrid between the Hernquist's TREECODE and the 
NCSA's ZEUS-3D hydrodynamics code. Thus, this program is very similar to the one used 
by Roettiger et al. (1996). The main differences are the use of a new Poisson solver for the 
creation of the gravitational potential, the inclusion of radiative cooling in the treatment of 
the ICM energy budget, the use of symmetry along the merger axis which allowed us to 
obtain better spatial resolution (20 kpc), and the ability to trace the gas in each subcluster 
with the help of passive scalars ( |G6mez et al. 1997b| ). 

These simulations start with two idealized clusters (main cluster and a subcluster) that 
are modeled as isothermal King spheres separated by ~ 4 Mpc that are allowed to interact 
under the effects of their own gravity. As the merger progresses, we observe the typical 
characteristics of merging clusters. For instance, we find that the gas is heated, compressed, 
and develops shocks. Moreover, the gas becomes clumpy and elongated as it tries to regain 
hydrostatic equilibrium with an evolving gravitational potential. All of these properties 
are similar to the ones identified by other numerical simulations of cluster mergers (e.g., 
Roettiger et al. 19961 , [Schindler fc Muller 1993|) . 



In order to compare these simulations with the actual data, we identified the simulation 
that most closely resembled the observed X-ray and optical morphology of the cluster. A 
merger between two subclusters with mass ratio of 4 to 1 about 1 Gyr after core crossing is 
our best matching simulation (compare Figure la and 2f with Figure 5). Figure 5 shows a 
map of the logarithm of the simulated gas density. We have also overlaid onto this Figure 
two contour maps which represent the surface density of the N-body particles (dashed 
contours). One of the contour plots shows the particles belonging to the main cluster (short 
dashes) while the other contour plot shows the distribution of the secondary cluster particles 
(long dashes). If we assume that the N-body particles represent the cluster galaxies, we 
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find very interesting morphological similarities. A close comparison between this simulation 
and the observations can explain several of the cluster properties. 

First, it will explain the apparent spatial, but not kinematical, coincidence of two 
distinct groups of galaxies. During this merger, especially if viewed close to the line of sight, 
one would expect that the galaxies of each pre-merger cluster would appear mixed in the 
plane of the sky. However, they would be separated in velocity space, as is the case for this 
cluster. 

Second, the X-ray substructure that we detect in the system can be explained by 
the elongations and dumpiness observed in the simulated gas density. These deviations 
from a relaxed, symmetric, and circular distribution are evidence of a system that is not 
in hydrostatic equilibrium. Evidence of this effect has been found in other samples of 
non- relaxed clusters (e.g., Mohr et al 1995). A similar result has been found in the analysis 
of synthetic X-ray images from numerical simulations of cluster mergers ( [Roettiger et al.| 

pi. 



Third, this simulation shows that gas belonging to the infalling subcluster is protected, 
and is not completely stripped by ram-pressure. The protective mechanism is the formation 
of a bow shock that surrounds the subcluster ( [Roettiger et al. 1996| ). This effect might 



explain the fact that the southern radio source is located within a region of X-ray emitting 
gas. In this model, the northern radio source is associated with the main cluster and hence 
to the main clump of X-ray emitting gas. 

Fourth, the coincidence between the position of the N-body particles with the position 
of the subcluster galaxies can explain the fact that the radio galaxies appear at rest with 
respect to their local subcluster. 

Fifth, we can explain the bending of the southern NAT from the interaction between 
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the radio jets with the dynamic ICM. In this model, the radio galaxy is located at rest with 
respect to its subgroup which is consistent with the optical data. However, the relative 
velocity of the ICM with respect to the galaxy could be large (~ 1000 km/s) as a result of 
a bulk flow of gas caused by the merger. 

Figure 6 shows a histogram of the relative velocity of N-body particles with respect to 
the gas velocity for the 4 to 1 merger simulation. For this plot, we have chosen the N-body 
particles located within the left clump of the secondary cluster in Figure 5. Next, if we 
assume that the N-body particles represent the cluster galaxies, Figure 6 shows the actual 
relative velocity of the subcluster galaxies surrounding the NAT with respect to the cluster 
ICM (again, compare Figure la with Figure 5). It is interesting to note that these velocities 
are distributed from ~ -1500 km/s to 2000 km/s. These are the velocities that are required 
to explain the jet bending (section 3.2.3). Therefore, despite the fact that the radio galaxy 
is at rest with respect to the other galaxies, the radio galaxy will not be at rest with respect 
to the local ICM. It is important to point out that we do not consider this simulation as a 
unique model for A578. We are only using these simulations as guidelines for the type of 
gas dynamics that are produced during cluster mergers. 



4.2. A1569 

At first glance, A1569 appears to be a bimodal X-ray cluster. The clumps of X-ray 
emission are associated with each of the two radio sources. Moreover, we have also identified 
two distinct groups of galaxies (N and S) also associated with the X-ray clumps and the 
radio galaxies. 

Are these two groups of galaxies gravitationally bound?. The velocity dispersion of 
the whole system (Sbi= 1618 km/s) appears to be too large for the X-ray luminosity and 
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for the overall optical richness of the cluster. We have applied a 2-body analysis to this 
system in order to determine if there is a projection angle that would find a bound solution 
for the system. In this analysis, we use the Newtonian criteria for gravitational binding 
( Peers et al. 1982]) and the binding ratio as described by Davis et al. (1995). A binding 
ratio < 0.39 will imply that there is a projection angle that will provide a bound solution 
for the system. In our case, we obtained a binding ratio of 2.06 which reveals that there 
is no bound solution for the system. It is important to note that we used mass estimates 
for the systems derived from the optical data and not from the X-ray observations during 
our analysis (Virial Mass for the northern group = 3.4 x 10 14 Mq, southern group = 1.1 x 
10 14 Mq). The system will be bound only if the mass estimates for the subgroups are gross 
underestimates of the real values (at least one order of magnitude) which seems unlikely. 

Once we divided the galaxies in these two groups, we concentrated our study on 
the southern group. Our analysis of the southern group of galaxies reveals important 
kinematical and X-ray properties. As is the case for A578, we consider that a similar merger 
model provides the best explanation for those properties. 

First, this model provides an explanation for the peculiar motion of the WAT. It 
has been suggested in the past ( Pegerle fc Hill 1994| , |Pinkney 1995Q that large peculiar 
motions observed for massive galaxies, especially for a WAT host, are the result of velocity 
contamination produced by a merger. 

Second, our X-ray morphological tests find evidence of substructure. Again, a look at 
Figure 5 can provide a reasonable explanation for the elongation and dumpiness detected 
in Al569's southern clump (Figure lb). This substructure suggests that the gas is not in 
hydrostatic equilibrium. 

Third, a merger model will also explain the similar morphologies between the spatial 
galaxy distribution and the X-ray emission (Figure 3f). Both distributions exhibit a E-W 
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elongation, which is a possible signature of a merger (see Figure 5). 

Fourth, we find that the WAT peculiar velocity is 215 ± 105 km/s. This velocity is 
not large enough to bend the WAT since ~ 1000 km/s is required according to our models 
(section 3.2.3). The observed WAT peculiar velocity is probably the result of velocity 
contamination produced by the galaxies from a merging subcluster. This merger could also 
explain the WAT bending as produced by the interaction between the dynamic ICM and 
the jet tails. This mechanism is similar to the one described for A578. Finally, Gomez et 
al. (1997) points out that there is a correlation between the local X-ray elongation and 
the direction of the bending of the jet tail. They explained that this alignment is another 
signature of the effects produced by a merger that bends and re-shapes the radio emitting 
plasma. 

5. SUMMARY AND CONCLUSIONS 

We obtained 74 new redshifts in A578 and find that 41 galaxies are cluster members. 
We also gathered 56 new velocity measurements in A1569 and confirmed that 54 of those 
galaxies are cluster members. Our analysis of these optical data, X-ray archival data, and 
radio data has revealed the following properties: 

First, both clusters show evidence of X-ray substructure. The significant substructure 
was detected as deviations from circular models of the cluster X-ray emission (clumps 
and elongations). The X-ray emission from A578 shows two elongations; one is in the 
N-S direction while the other is in the E-W direction (Figure la). A1569 appears to be a 
bimodal X-ray cluster with an elongated (E-W direction) southern clump (Figure lb). 

Second, our 1-D, 2-D, and 3-D statistical tests find optical substructure in both 
clusters. The KMM algorithm successfully identifies two subgroups of galaxies present in 
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these sytems. Thus, both clusters are divided into two subgroups each. Each subgroup 
contains one extended radio source. 

Third, the positions of the galaxy subgroups coincide with the position of the local 
X-ray clumps. This seems to suggests that these subgroups are composed of galaxies and 
gas. A more detailed analysis reveals that the two subgroups in A578 are bound and 
coincide spatially but not kinematically. On the other hand, our 2-body analysis of A1569 
finds that its two subgroups are not gravitationally bound. 

Fourth, it is interesting that the southern NAT in A578 appears to be at rest with 
respect to its subgroup galaxies. However, our model of ram induced jet bending requires a 
relative velocity of the jet with respect to the ICM of ^ 1000 km/s in order to reproduce 
the observed morphology. So, what causes the NAT bending? 

Fifth, the WAT in A1569 exhibits a small peculiar velocity (~ 200 km/s) which is not 
enough to explain the bending of the WAT according to a ram-pressure model. 

We believe that a cluster-subcluster merger model provides a reasonable explanation 
for the observed optical and X-ray properties of these clusters and for the bending of their 
radio sources. During a merger, the cluster gas will be disrupted and shocked as it tries to 
regain hydrostatic equilibrium. Morphologically, the gas will show clumps and elongations 
that are similar to the ones seen in our observations of these clusters. This dynamic 
environment also produces bulk flows of gas 1000 km/s) that in turn are responsible for 
the bending of the radio tails. Thus, we believe that a merger is the best mechanism to 
explain the presence of extended and bent radio sources in these clusters. Therefore, we 
have confirmed our suspicions that extended radio sources act as beacons of merger activity 
in these two clusters. 

Finally, we would like to point out that these two clusters were chosen simply because 
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they contained galaxies with extended radio emission and had available ROSAT PSPC 
data. We are in the process of performing a similar analysis on a larger sample (~ 15) of 
galaxy clusters that contain extended radio sources and on a control sample of radio-quiet 
galaxy clusters. In this way, we will be able to use a larger sample for determining if 
extended radio sources are more likely to be found in clusters that are undergoing mergers. 
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Figure Captions 

Fig la: Overlay of a 20 cm radio map grey scale onto ROSAT X-ray surface brightness 
contours. The contour map shows the adaptive smoothed X-ray emission from A578. The 
contour levels are 0.06, 0.10, 0.20, 0.30, 0.50, 0.60, 0.70, 0.80, and 1.00 counts arcmin- 2 . 

Fig lb: As in Figure la) for the cluster A1569. The contour levels are 0.06, 0.12, 0.24, 
0.36, 0.48, 0.60, 0.84, 1.08, and 1.32 counts arcmin^ 2 . 

Fig 2a: Velocity histogram of all the galaxies observed in the A578 field. The binsize 
is 200 km/s. We have also marked (crosses) the positions of the radio galaxies. 

Fig 2b: 2-D velocity-coded plot of all the galaxies observed in the A578 field. The 
triangles correspond to the galaxies with velocities between 10,000 and 24,000 km/s. The 
squares correspond to the 24,000-28,000 km/s velocity range. Finally, circles correspond 
to the galaxies with velocities > 28,000 km/s. The stars mark the positions of the radio 
galaxies. 

Fig 2c: Velocity histogram of the 41 possible galaxy members for A578. We have also 
marked (crosses) the positions of the radio galaxies. 

Fig 2d: 2-D velocity-coded plot of the 41 possible galaxy members for A578. The 
circles represent the galaxies with velocities between 25,000 and 27,000 km/s. The squares 
show the positions of the galaxies with velocities between 27,000 and 29,000 km/s. Again, 
the stars mark the position of the radio galaxies. 

Fig 2e: 2-D plot of the result of the KMM objective partition of the galaxies. We 
have also marked the position of the radio galaxies (stars). The squares correspond to the 
members of the North group whereas circles represent the galaxies belonging to the South 
group. 
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Fig 2f: Overlay of an X-ray map grey scale onto galaxy surface density contours for 
A578. The continuous contours represent the surface density of the northern subgroup with 
the following levels: 0.014, 0.028, 0.048, 0.067, 0.084, 0.106, and 0.12 galaxies arcmin- 2 . 
The dashed contours correspond to the surface density of the southern galaxy subgroup. 
Their contour levels are 0.01, 0.02, 0.03, 0.04, 0.06, 0.09, 0.1, 0.18, 0.23, 0.29, 0.34, 0.4, 
and 0.5 galaxies arcmin" 2 . Note how the elongation of the southern subgroup mimic the 
elongation of the central X-ray emission. 

Fig 3a: Velocity histogram of all the galaxies observed in the A1569 field. The binsize 
is 200 km/s. We have also marked (crosses) the positions of the radio galaxies. 

Fig 3b: Velocity histogram of the 41 possible galaxy members for A1569. We have 
also marked (crosses) the positions of the radio galaxies. 

Fig 3c: 2-D velocity-coded plot of all the galaxies observed in the A1569 field. The 
squares show the positions of the galaxies with velocities between 19,000 and 22,000 km/s. 
The circles represent the galaxies with velocities between 22,000 and 25,000 km/s. The 
stars mark the position of the radio galaxies. 

Fig 3d: 2-D plot of the result of the KMM objective partition of the galaxies. We 
have also marked the position of the radio galaxies (stars). The squares correspond to 
the members of the southern group whereas circles represent the galaxies belonging to the 
northern group. 

Fig 3e: Velocity histogram of the result of the KMM objective partition of the galaxies 
in A1569. The binsize is 200 km/s. We have also marked the positions of the radio galaxies 
(squares). The dashed lines represent the velocity distribution of the galaxies that belong 
to the southern subgroup. 

Fig 3f: Overlay of an X-ray map grey scale onto a galaxy surface density contour map. 
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The contour map represents the surface density of all the galaxies within 19,000-25,000 
km/s. The contour levels are 0.15, 0.25, 0.4, 0.55, 0.7, and 0.85 galaxies arcmin -2 . Note 
how there are two main groups of galaxies that coincide with the two main peaks of X-ray 
emission. It also appears that these two groups are connected by a filament of X-ray 
emitting gas and galaxies. 

Fig 4: Relative galaxy velocity and number density as a function of (3j for the NAT in 
A578 and the WAT in A1569. The thick lines correspond to the density while the thin lines 
correspond to the galaxy velocity. The efficiency of the conversion between the jet energy 
and the radio luminosity is indicated by the percentages. 

Fig 5: Overlay of three contour plots. The continuous contour plot represents the gas 
density from a 4 to 1 mass ratio merger simulation 1 Gyr after core crossing. These contours 
have the following levels: 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, and 19 
times 0.00023 cm -3 . The long dashed contours represent the surface density of N-body 
particles that belong to the main cluster with the following contour levels with respect to 
the peak value: 100%, 84%, 42%, 21%, and 11%. This main cluster had a pre-merger mass 
4 times greater than the secondary cluster. The short dashed contours represent the surface 
density of the particles belonging to the secondary cluster with the following contour levels 
in units of the peak value: 100%, 77%, 52%, and 26%. Note how the spatial elongations of 
these two subgroups follow the direction of the merger (from left to right) and how they 
mimic the distribution of the galaxies in A578. Each tic mark is spaced by 20 kpc. 

Fig 6: Histogram of the relative velocity of the gas with respect to the N-body 
particles for a 4 to 1 merger simulation. The particles correspond to one of the subgroups 
of the secondary cluster (the rightmost subclump of particles). The time corresponds to 1 
Gyr after core crossing. If these particles are similar to the galaxies of the cluster, then this 
plot shows the effective relative velocity of the gas with respect to the galaxies that can be 
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responsible for the bending of the radio jets (NAT and WAT). 
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TABLE 1. Observation Log 



Cluster Date of Observation 



A578 29-31 January 1995 
19-20 October 1995 
17-19 February 1996 
5-7 October 1996 

A1569 17-19 February 1996 
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TABLE 2 A. Galaxies in the A578 field 



Galaxy ID 


KA(zOOO) 


UEO(2000J 


Heliocentric Velocity 


error 


1UK value 


1 


7 24 54.09 


66 59 09.90 


26150 


35 


12.8 


2 


7 24 56.26 


66 59 09.10 


26106 


51 


8.6 


3 


7 25 08.93 


66 59 24.80 


25842 


33 


13.8 


4 


7 25 23.11 


66 58 31.90 


26456 


51 


8.5 


5 


7 25 01.99 


66 56 32.10 


27709 


45 


9.8 


6 


7 24 03.48 


66 59 25.20 


25797 


89 


4.5 


7 


7 24 38.96 


67 00 20.00 


26057 


94 


4.2 


8 


7 24 32.34 


67 01 22.10 


26106 


89 


4.5 


9 


7 24 33.99 


67 03 39.00 


27387 


86 


4.7 


10 


7 25 12.10 


67 04 26.50 


27613 


52 


8.4 


11 


7 24 58.82 


67 02 32.30 


28070 


83 


4.9 


12 


7 25 13.29 


67 02 19.10 


27578 


90 


4. 


13 


7 25 14.77 


67 00 47.70 


39272 


56 


7.8 


14 


7 25 44.81 


66 58 25.10 


11883 


42 


10.5 


15 


7 25 42.31 


66 54 50.40 


26938 


69 


6.1 


16 


7 24 55.22 


66 53 33.10 


26658 


40 


11.1 


17 


7 23 52.63 


66 55 38.50 


26604 


103 


3.7 


18 


7 26 47.88 


67 02 12.90 


27415 


61 


6.9 


19 


7 27 17.22 


66 59 03.90 


27389 


85 


4.7 


20 


7 27 23.87 


67 02 19.40 


26334 


53 


8.2 


21 


7 27 34.63 


67 03 49.30 


26381 


77 


5.3 


22 


7 27 27.43 


67 08 50.70 


22407 


52 


8.4 


23 


7 27 25.18 


67 12 14.30 


40088 


86 


4.7 


24 


7 25 28.27 


67 18 07.90 


39702 


82 


4.9 


25 


7 24 26.34 


67 22 15.60 


19063 


85 


4.7 


26 


7 24 05.44 


67 14 43.30 


23735 


80 


5.1 


27 


7 24 19.20 


67 06 58.80 


27002 


69 


6.1 


28 


7 23 55.52 


67 07 37.30 


42110 


85 


4.7 


29 


7 23 03.95 


67 11 35.10 


42162 


72 


5.8 


30 


7 22 44.92 


67 17 49.10 


28488 


77 


5.4 


31 


7 22 35.31 


67 12 54.90 


42393 


80 


5.1 


32 


7 22 09.55 


67 11 12.60 


23759 


95 


4.1 


33 


7 22 45.08 


67 01 20.60 


40600 


71 


5.8 


34 


7 24 15.76 


66 49 06.30 


26094 


65 


6.5 


35 


7 23 23.80 


66 50 41.70 


39299 


84 


4.8 


36 


7 23 29.82 


66 49 36.50 


39406 


95 


4.1 


37 


7 21 30.47 


66 50 44.60 


40651 


79 


5.1 


38 


7 22 01.33 


66 37 37.70 


21023 


72 


5.7 


39 


7 27 45.46 


66 52 52.70 


26825 


81 


5.0 


40 


7 26 37.99 


66 53 14.10 


26167 


83 


4.8 


41 


7 26 50.31 


66 54 17.30 


27685 


72 


5.8 


42 


7 25 04.47 


66 59 29.20 


25803 


101 


3.8 


43 


7 24 51.92 


66 59 23.60 


26431 


80 


5.1 


44 


7 24 54.27 


66 58 31.60 


28134 


68 


6.1 


45 


7 24 38.10 


66 58 11.90 


28017 


75 


5.5 


46 


7 25 14.51 


66 57 56.90 


16449 


104 


3.7 


47 


7 25 14.93 


66 57 41.40 


25252 


103 


3.7 


48 


7 25 38.78 


66 59 31.70 


33957 


102 


3.8 



TABLE 2 A. (continued) 



Galaxy ID 


RA(2000) 


DEC(2000) 


Heliocentric Velocity 


error 


TDR value 


49 


7 25 40 46 


66 59 26 00 


23714 


101 


3.8 


50 


7 26 08 51 


66 58 47 90 


19663 


108 


3.5 


51 


7 25 1 3 43 


66 55 49 60 


28213 


86 


4.6 


52 


7 25 00 75 


66 56 23 20 


26386 


78 


5.2 


53 


7 24 25 22 


66 55 1 80 

\J \J 'J 'J 1W.UW 


26344 


90 


4.4 


54 


7 23 54 28 


66 55 45 30 


11412 


105 


3.6 


55 


7 24 18.84 


67 01 49.90 


26494 


82 


5.0 


56 


7 24 12.21 


67 02 24.50 


40636 


100 


3.9 


57 


7 24 29.32 


67 02 18.60 


39071 


97 


4.0 


58 


7 24 37.13 


67 03 10.30 


39184 


99 


3.9 


59 


7 24 32.97 


67 05 34.40 


27466 


97 


4.0 


60 


7 24 59.68 


67 03 29.20 


8421 


103 


3.7 


61 


7 26 03.90 


67 02 29.60 


26076 


77 


5.4 


62 


7 25 51.58 


67 07 13.90 


41884 


95 


4.1 


63 


7 25 41.50 


67 08 43.20 


28055 


63 


6.7 


64 


7 25 51.81 


67 17 27.90 


25800 


92 


4.3 


65 


7 21 18.64 


67 11 25.90 


20694 


89 


4.4 


66 


7 21 51.08 


67 04 32.30 


3551 


97 


4.0 


67 


7 21 41.59 


67 02 33.70 


14785 


90 


4.4 


68 


7 22 03.05 


66 59 28.50 


39539 


99 


3.9 


69 


7 21 55.29 


66 59 24.30 


20873 


93 


4.2 


70 


7 22 02.42 


66 55 14.40 


39013 


67 


6.2 


71 


7 23 14.23 


66 55 05.70 


26578 


84 


4.8 


72 


7 22 12.07 


66 51 17.60 


39270 


96 


4.1 


73 


7 23 12.30 


66 45 16.00 


26941 


95 


4.1 


74 


7 23 24.75 


66 45 02.10 


26781 


88 


4.5 



TABLE 2B. Galaxies in the A1569 field 



Galaxy ID 


KA(zOOO) 


UEO(2000J 


Heliocentric Velocity 


error 


1UK value 


1 


12 36 08.15 


16 38 28.7 


23698 


73 


5.7 


2 


12 36 25.80 


16 32 18.4 


20525 


61 


7.0 


3 


12 36 14.69 


16 32 55.7 


20854 


57 


7.6 


4 


12 36 37.49 


16 33 44.9 


22111 


73 


5.7 


5 


12 36 37.29 


16 38 57.6 


20832 


73 


5.7 


6 


12 35 59.74 


16 35 33.6 


21191 


90 


4.4 


7 


12 36 37.98 


16 28 43.7 


31526 


62 


6.8 


8 


12 36 35.66 


16 26 36.8 


23611 


88 


4.5 


9 


12 35 25.36 


16 29 59.6 


23809 


75 


5.5 


10 


12 36 17.72 


16 24 26.7 


21280 


50 


8.7 


11 


12 36 08.52 


16 38 50.2 


22711 


38 


11.9 


12 


12 36 07.95 


16 38 51.5 


24183 


39 


11.6 


13 


12 36 26.33 


16 36 58.3 


24450 


75 


5.5 


14 


12 36 21.50 


16 32 49.3 


20300 


70 


5.9 


15 


12 36 19.12 


16 32 46.5 


20149 


32 


14.4 


16 


12 36 24.53 


16 32 22.7 


19971 


62 


6.9 


17 


12 36 27.85 


16 32 04.5 


20361 


85 


4.7 


18 


12 36 13.78 


16 32 42.2 


20824 


53 


8.2 


19 


12 36 05.12 


16 32 19.3 


21145 


67 


6.3 


20 


12 36 32.22 


16 39 23.1 


23539 


82 


4.9 


21 


12 35 54.28 


16 39 39.0 


24266 


106 


3.6 


22 


12 35 45.57 


16 35 07.6 


31652 


99 


3.9 


23 


12 35 39.79 


16 26 26.6 


20704 


44 


10.1 


24 


12 36 22.43 


16 29 47.8 


20852 


64 


6.6 


25 


12 36 30.72 


16 31 19.2 


19469 


75 


5.5 


26 


12 35 24.28 


16 35 11.4 


21186 


68 


6.2 


27 


12 37 12.84 


16 34 33.7 


24280 


90 


4.4 


28 


12 36 10.64 


16 38 29.2 


22604 


76 


5.4 


29 


12 36 05.73 


16 38 28.8 


23550 


51 


8.5 


30 


12 36 15.82 


16 35 35.3 


20759 


69 


6.0 


31 


12 36 20.84 


16 38 36.5 


24166 


72 


5.8 


32 


12 36 24.78 


16 35 22.7 


24329 


62 


6.9 


33 


12 36 27.05 


16 33 47.4 


20359 


57 


7.6 


34 


12 36 27.54 


16 32 19.3 


20707 


77 


5.3 


35 


12 36 11.39 


16 44 44.6 


23511 


84 


4.8 


36 


12 36 14.67 


16 39 06.0 


23078 


61 


7.0 


37 


12 35 42.60 


16 38 58.5 


23563 


64 


6.6 


38 


12 35 50.13 


16 33 15.6 


22749 


72 


5.8 


39 


12 35 47.02 


16 27 51.1 


20940 


79 


5.2 


40 


12 35 51.91 


16 25 51.4 


23525 


97 


4.0 


41 


12 36 18.61 


16 28 57.0 


20677 


84 


4.8 


42 


12 36 26.72 


16 31 12.1 


20553 


67 


6.3 


43 


12 36 26.72 


16 31 43.5 


20117 


37 


12.0 


44 


12 36 24.84 


16 31 59.3 


20823 


88 


4.5 


45 


12 36 44.14 


16 32 02.6 


24058 


75 


5.5 


46 


12 36 42.27 


16 26 20.8 


20057 


84 


4.8 


47 


12 36 26.04 


16 36 31.7 


24649 


103 


3.7 


48 


12 36 36.15 


16 34 16.0 


21075 


92 


4.3 



TABLE 2B. (continued) 



Galaxy ID 


RA(2000) 


DEC(2000) 


Heliocentric Velocity 


error 


TDR value 


49 


12 36 04.75 


16 37 32.9 


24852 


104 


3.7 


50 


12 35 44.49 


16 26 31.6 


20816 


74 


5.6 


51 


12 36 17.33 


16 22 38.4 


21541 


62 


6.8 


52 


12 36 24.30 


16 23 20.9 


23725 


47 


9.3 


53 


12 36 37.60 


16 36 30.0 


21359 


65 


6.5 


54 


12 36 03.68 


16 44 59.0 


24005 


80 


5.1 


55 


12 35 43.12 


16 44 19.0 


24035 


85 


4.7 


56 


12 35 36.37 


16 40 23.7 


20906 


87 


4.6 
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TABLE 3 A. Spectral Fitting Regions 



Cluster 


RA(1950) 


DEC(1950) 


Source 


Background 








Radius a 


Radius/ Width a 


A578 


07 19 55 


67 05 20 


5.7 


21.75/2 


A1569 


12 33 57 


16 48 39 


4.9 


12.5/3 



a The units are arcminutes. 



TABLE 3B. Spectral Fitting Results' 



Cluster 




Net Counts/ 
Error 


Hi d 
- 1 - spectral 

(keV) 


xil 

d. o. f. 


' 1 d 

i-C-stat 

(keV) 


A578 


20.64 


201/16 


1 7+ 1 ' 5 

1 -' -0.5 


0.88/6 


1 7+ 1 ' 7 

1 -' -OA 


A1569 


20.41 


210/16 


i 4+° 9 


0.8/4 


i 4+ 1 



a All the fits were performed with a cosmic abundance of 0.3. 
b The values for log(N#) are from Stark et al. (1992). 
c Net counts is the number of background-subtracted counts. The 
error in the counts corresponds to the la value. 

d The errors are 68% confidence level for 1 free parameter. 



TABLE 4. Optical Data for A578 and A1569 



cluster 


CWtkms" 1 ) 


Ss/(kms - 1 ) 


~N galaxies 


A578 
A1569 


26,772t^° 
22,058t*j 9 


793±?f 
1,618±H 


41 
54 



TABLE 5. Statistical Analysis for A578 



test 


value si f 


J"n ifi ran cp 


A 


0.85 




U 


3.97 


0.35 


W 


0.95 


0.07 


Bl 


0.34 


0.16 


B2 


2.08 


0.05 


J3 1-J3Z 


o.ol 


U. 10 


i 


0.96 


0.15 


DIP 


0.05 


0.5 


KS 


0.77 


0.25 


KSv 


1.43 


0.10 


W 2 


0.14 


0.03 


u 2 


0.13 


0.03 


A 2 


0.81 


0.04 


TI 


0.86 


0.15 


AI 


0.86 


0.05 


normal 




0.12 


skewness 




0.16 


kurtosis 




0.09 


2D tests 


Lee 2D 




0.30 


FEL 




0.50 






0.26 


AST 




0.62 


3D tests 


Lee 3D 




0.34 


A 




0.57 


a 




0.02 


a var 




0.25 


e var 




0.11 



Notes to Table 5. 

We are only listing the significance 
for the 2D and 3D tests evaluated 
from 500 Monte Carlo simulations 
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■3000 -2000 -1000 1000 2000 
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TABLE 6 A. A578 Optical KMM Subgroups Data 



parameter Value for the Northern Subgroup Value for the Southern Subgroup 

RA(2000) 7 25 40.39 7 24 42.66 

DEC(2000) 67 01 18.6 66 57 52.8 

~N galaxies 18 23 

Cfi/fkms- 1 ) 27,5791^7 26,266tg3 

S B/ (kms- 1 ) 565±H1 385±|? 



TABLE 6B. KMM Galaxy Assignments 

Group Galaxies 

N 1, 2, 3, 4, 6, 7, 8, 16, 17, 27, 34, 40, 42, 43, 47, 52, 53, 55, 61, 64, 71, 73, 74 
S 5, 9, 10, 11, 12, 15, 18, 19, 20, 21, 30, 39, 41, 44, 45, 51, 59, 63 



TABLE 7. Statistical Analysis for A1569 



test 


value sij 


jjnificance 


A 


0.93 




U 


3.32 


0.001 


W 


0.87 


0.001 


Bl 


0.22 


0.24 


B2 


1.46 


0.00 


B1-B2 


63.44 


0.00 


I 


1.07 


0.15 


DIP 


0.09 


0.99 


KS 


1.45 


0.01 


KSv 


2.82 


0.01 


W 2 


0.50 


0.001 


u 2 


0.50 


0.001 


A 2 


2.78 


0.001 


normal 




0.03 


skewness 




0.13 


kurtosis 




0.001 


TI 


0.72 




AI 


0.52 




2D tests 


Lee 2D 




0.03 


FEL 




0.14 






0.002 


AST 




0.42 


3D tests 


Lee 3D 




0.04 


A 




0.14 


a 




0.19 


a var 




0.44 


e var 




0.01 



Notes to Table 7. 

We are only listing the significance 
for the 2D and 3D tests evaluated 
from 500 Monte Carlo simulations 



TABLE 8 A. A1569 Optical KMM Subgroups Data 



parameter Value for the Northern Subgroup Value for the Southern Subgroup 

RA(2000) 12 36 13.07 12 36 15.65 

DEC(2000) 16 36 21.09 16 31 48.27 

~N galaxies 25 29 

Cfi/fkms- 1 ) 23,7821^5 20,740^? 

S B/ (kms- 1 ) 622±ii£ 433tf 7 



TABLE 8B. KMM Galaxy Assignments 



Group Galaxies 



N-E 1, 4, 8, 9, 11, 12, 13, 20, 21, 27, 28, 29, 31, 32, 35, 36, 37, 38, 40, 45, 47, 49, 52, 54, 55 

S-W 2, 3, 5, 6, 10, 14, 15, 16, 17, 18, 19, 23, 24, 25, 26, 30, 33, 34, 39, 41, 42, 43, 44, 46, 48, 50, 51, 53, 56 





Table 9 


K X-ray 


Substructure 


Analysis 


Cluster 


a 


8 


Probability 


Notes 




(1950) 


(1950) 






A578 


07 20 15 


67 05 00 


0.05 


Western elongation 




07 19 35 


67 04 30 


0.06 


Eastern elongation 




07 20 02 


67 02 30 


<0.01% 


S-W elongation 


A1569 


12 33 47 


16 49 30 


0.003 






12 34 07 


16 50 30 


<0.01% 






12 34 00 


16 47 00 


0.007 


S-E diffuse structure 



TABLE 10. Pressure Analysis 



Cluster 


a 

(1950) 


8 

(1950) 


Minimum Radio 
Source Pressure 
(10~ 13 dyne/cm 2 ) 


Deprojection 
Thermal Pressure 
(10~ 13 dyne/cm 2 ) 


Projected Distance 
from X-ray Center 
(kpc) 


(0.5 to 2.05 keV) 
(10 43 erg/s) 


A578 


07 19 49.0 


67 05 01 


55 


20±2 


50 


1.1 ± 0.1 




07 20 00.0 


67 02 48 


18 


10±2 


200 






07 20 05.0 


67 02 19 


10 


9±2 


250 




A1569 


12 33 59.5 


16 47 29 


3 


46±10 


85 


1.7±0.3 



